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Abstract. We present a new mode of hydrogen burning on neutron stars (NSs) 
called diffusive nuclear burning (DNB). In DNB, the burning occurs in the ex¬ 
ponentially suppressed tail of hydrogen that extends to the hotter regions of the 
envelope where protons are readily captured. Diffusive nuclear burning changes 
the compositional structure of the envelope on timescales ~ 10^“^ yrs, much 
shorter than otherwise expected. This mechanism is applicable to the physics of 
young pulsars, millisecond radio pulsars (MSPs) and quiescent low mass X-ray 
binaries (LMXBs). 


1. Introduction 

The composition of NS envelopes affects their cooling (Potekhin, Chabrier & Yakovlev 
1997) and thermal radiation (Romani 1987). The amount of material needed to change 
the spectral profile of the thermal radiation from the surface of a NS is miniscule 
(~ IO^^^Mq). Such a small amount of contamination could easily be produced from 
spallation (Bildsten, Salpeter & Wasserman 1992) of fallback material (Woosley & 
Weaver 1995) soon after a supernova explosion. 

Recent X-ray observations of the thermal spectrum from young radio pulsars have 
yielded tantalizing clues of their photospheric makeup (see Pavlov, Zavlin & Sanwal 
2002 for a review). The thermal emission from young NSs can be fit with two models: 
magnetic hydrogen atmospheres or blackbody atmospheres. Both atmospheres fit the 
thermal spectra equally well, however they yield different effective temperatures and 
different solid angles for the emission area. The model which yields a more reasonable 
radius for the NS at the preferred distance is then taken as the favored model (Pavlov 
et al. 2002). For radio pulsars younger than ^ 10'^“®yrs (e.g. Vela), the radius is 
reasonable for a magnetic hydrogen atmosphere model (Pavlov et al. 2001). For pulsars 
older than ~ 10^“^ yrs (e.g. PSR B0656-I-14), a blackbody model is favored (see Pavlov 
et al. 2002 and references therein). This suggests a possible evolution of hydrogen to 
more blackbody like elements on a timescale of 10^“^ yrs, and these indications have 
motivated our work. 

At the photosphere, the local temperature (T ~ 10® K) and density [p ~ lgcm“®) 
are too low to allow any significant nuclear evolution over 10^“® yrs. However at a depth 
1 m underneath the photosphere, the temperature is roughly two orders of magnitude 
greater (T ~ 10® K). Consider a NS envelope of hydrogen on carbon (we use C here as 
our fiducial proton capturing nucleus, results for other nuclei are in Chang & Bildsten 
2003; CB03 hereafter) as shown in Figure 1. In diffusive equilibrium, the separation 
between the hydrogen and carbon is not strict. Rather a diffusive tail of hydrogen pene¬ 
trates deep into the carbon layer. Protons easily reach a depth where the temperature is 
sufficiently high for rapid capture onto C. Since the hydrogen’s diffusive tail is exponen- 
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Figure 1. Schematic of a H/C envelope in diffusive equilibrium. The diffusive 
tail of hydrogen extends deep into the carbon, reaching temperatures where 
the hydrogen rapidly captures onto carbon, forcing depletion of the hydrogen 
layer as protons diffuse down to the burning layer. 


tially suppressed, there will be a region where proton captures are peaked, which we call 
the burning layer. The consumption of the hydrogen by burning will drive the diffusive 
tail out of equilibrium, which sets up a diffusive current of hydrogen that flows down 
from the hydrogen layer into the burning layer. We refer to this burning mechanism 
(first mentioned by Chui & Salpeter 1964 and initially calculated by Rosen 1968) as 
diffusive nuclear burning (DNB; CB03). Over time, given that there is no interstellar 
accretion or other processes to refresh the hydrogen layer, the H is depleted. 

This scenario presents several competing timescales. The first is the proton capture 
timescale in the burning layer, Tnuc- The second is the time it takes protons to diffuse 
into the burning layer, Tdig. In the case where Tdis ^ t^ucj the diffusive tail is always 
in equilibrium. In the opposite limit, the diffusive tail is modified by the burning. For 
simplicity we will only discuss the first case here. 

The equilibrium structure of the diffusive tail is calculated from hydrostatic balance 
for each ion, 


= -m {Aimpg - Z^eE ), (1) 

= -Tie {meg + eE ), (2) 

where Pi, rii, Ai, Zi are the pressure, number density, atomic number and charge of the 
Fth ion species and E is the upward pointing electric field found by demanding charge 
neutrality, Ue = 'Y^ntZi. The thermal structure is determined from the constant flux 
equation. 


dr 
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where k is the opacity and Te is the effective temperature. Given appropriate micro¬ 
physics, the equilibrium thermal and compositional structure can be calculated (see 
CB03 for additional details). The microphysics we have chosen in this paper are valid 
for the non-magnetic case, B < 10® G. 
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Figure 2. Differential hydrogen column burning rate taking into account p- 
p capture and p + C capture. The bottom graph shows the number fraction 
(solid line) and temperature (dotted line). This model has yn = 100gcm“^ 
and Te = 8 X 10^ K. The integrated burning rate for this model is yn/ tco\ = 
0.24gcm“2yr“\ giving TcoI = 417yrs. 


The local hydrogen burning rate is mpUH/rnuc, hence the total hydrogen burning 
rate per area on a NS, C//, is 


^ _ f nurup 
Tcol J TH{nH,nc,T) 


( 4 ) 


where yn = J mpUndz is the integrated column of hydrogen and Tcoi is the characteristic 
burning time for that column. In Figure 2, we show the equilibrium structure for a 
fiducial NS model with Tg = 8 x 10^ K and the total burning rate per log column, 
y = P/g. The Gaussian peak in the burning rate traces out the burning zone which is 
centered around a column of y « 10®gcm“^, which is about 20 cm below where most 
of the H resides. 

We relate the burning rate per column into a total mass burning rate Mdnb = 
47ri?^Ci/. In Figure 3, we plot the characteristic burning time Tcoi and associated Mdnb 
as a function of the total hydrogen column, yn- For a given envelope with an initial 
column of hydrogen at a fixed core temperature, the evolution follows the curve up to the 
photosphere at yu ~ lgcm“^. Because the evolution follows a simple power law, the 
lifetime for any given envelope is completely determined by the characteristic burning 
time. Tool, at the photosphere. The curve for a central temperature of Tg = 5 x 10^ K is 
cut off at large columns since the energy release from proton captures is comparable to 
the flux leaving the envelope. Hence our assumption of constant flux is violated. What 
is also notable about Figure 3 is that the curves follow a power law dependence which 
we derive in CB03. 

Future directions for DNB are including the microphysics for magnetic NSs, calcu¬ 
lating DNB without the constraint Tdiff ^ Tnuc and allowing for an intervening helium 
layer. DNB allows for young neutron stars to have atmospheres other than hydrogen 
soon after birth. For instance, NSs with carbon, nitrogen or oxygen photospheres can 
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Figure 3. Characteristic burning time, TcoI, and total mass burning rate, 
-Mdnb for a NS radius of 10 km, as a function of yn for different fixed core 
temperatures and an underlying carbon layer. For each model, we list the 
logarithmic core temperature and associated logarithmic effective temperature. 


easily exist after burning off any initial H. DNB may help explain the recent observation 
of magnetic oxygen lines on 1E1207+56 (see De Luca’s review in this volume, also see 
Hailey & Mori 2002). 
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